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We consider the signatures of a domain wall produced in the spontaneous symmetry breaking 
involving a dilaton-like scalar field coupled to electromagnetism. Domains on either side of the wall 
exhibit slight differences in their respective values of the fine-structure constant, a. If such a wall 
is present within our Hubble volume, absorption spectra at large redshifts may or may not provide 
a variation in a relative to the terrestrial value, depending on our relative position with respect to 
the wall. This wall could resolve the "contradiction" between claims of a variation of a based on 
Keck/Hires data and of the constancy of a based on VLT data. We derive the properties of the 
wall and the parameters of the underlying microscopic model required to reproduce the possible 
spatial variation of a. We discuss the constraints on the existence of the low-energy domain wall 
and describe its observational implications concerning the variation of the fundamental constants. 
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There are very few observations which can be directly 
and unambiguously related to new physics. The study of 
relative wavelength shifts in quasar absorption spectra at 
high redshift is indeed one of them as systematic achro- 
matic shifts in these spectra can be attributed to changes 
in fundamental constants, and in particular in the fine- 
structure constant, a. This would most certainly call for 
physics beyond the standard model. Study of the varia- 
tion of constants on cosmological scales is also the best 
way to test the equivalence principle on cosmological and 
astrophysical scales [l| . It opens a window on deviations 
from general relativity on scales where it is necessary to 
introduce dark energy and dark matter and on which we 
have very little constraint on the validity of general rela- 
tivity 0. 

Claims of a variation in a from observations of quasar 
absorption spectra using the many multiplet method Q 
had sparked an enormous amount of theoretical activity 
in attempts to explain a temporal variation in the fine 
structure constant 0411 ■ ^ confirmed, the Keck/Hires 
data which yielded a statistically significant trend indi- 
cating Aa/a = (— 0.54±0. 12) x 10~ 5 over a redshift range 
0.5 < z < 3.0 (the minus sign indicates a smaller value of 
a in the past) could indeed point to new physics. How- 
ever, subsequent studies based on VLT data using the 
same method have shown Aa to be consistent with zero 
9, 10]. These results remain somewhat controversial [llj ]. 



If the low energy constants of physics depend on some 
dynamical scalar field, (f>, they become dynamical and 
may well be space-time dependent. On cosmological 
scales, it is usually thought that the time variation dom- 
inates over spatial fluctuations, as suggested by most 
models. The reasoning here is straightforward. For a 
scalar field coupled to electromagnetism, the Lagrangian 
contains a term F^ V F ^ , F^ being the Faraday 

tensor and Bp an arbitrary function of cj>. This will nec- 
essarily induce a coupling to matter which is generated 
radiatively if not present at the tree level (see below). 
The equation of motion for the scalar field simply takes 
the form 



dV ci 



0. 



(1) 



where V e s includes the self interactions of </> as well as 
any couplings to matter, so that it may depend on the 
local energy density of matter. For example, should the 
Lagrangian contain a term B^(cj))mNNN, then the cou- 
pling to matter is effectively density dependent, which 
could serve as the source of spatial variations through 



U(f> + mfa = B' N (cj})p N , 



(2) 



where is the scalar mass and pat is the baryon en- 
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ergy density. 1 However, the density dependent shifts 
from the homogeneous solution are typically extremely 
small except perhaps in the vicinity of a neutron star [13l — 
[l5j . In contrast, temporal variations are relatively easy 
to achieve particularly over cosmological time scales, as 
long as the field remains light. 

However, there have been a series of recent puz- 
zling observational results. First, the combined positive 
Keck/Hires and negative VLT results for a change in a 
could be interpreted as a dipole in the spatial distribu- 
tion of a [16l - ll8j . Then, it has also been claimed that the 
ratio rrhn/rrie has a small spatial variation in the Milky 
Way . While caution should still rule the day (the 
positive result for a variation in a has yet to be con- 
firmed independently and may still be due to systematic 
effects 20] ) , it is an intriguing possibility with potentially 
very interesting interpretations. 

As we have argued above, spatial variations are ex- 
pected to be much smaller than a time variation. In- 
deed, if the field that triggers the variation of the con- 
stant is light during inflation, it would have developed 
super-Hubble fluctuations of quantum origin, with an al- 
most scale invariant power spectrum. The constants de- 
pending on such a field must also fluctuate on cosmolog- 
ical scales and have a non- vanishing correlation function. 
This possibility is however constrained [2l|, and would 
not be dipole in nature. 

Another possibility would be that the Copernican prin- 
ciple is not fully satisfied. Then, the background value of 
4> would depend e.g. on r and t for a spherically symmet- 
ric spacetime (such as a Lemaitre-Tolman-Bondi space- 
time) . This could give rise to a dipolar modulation of the 
constants if the observer (us) is not located at the center 
of the universe. Such a cosmological dipole should also 
reflect itself on other cosmological observations such as 
the cosmic microwave background (CMB) anisotropics, 
which does not seem to match with the required dipole 
in a [TH, [TtJ • Note also that such large scale deviations 
from homogeneity are constrained observationally [22j. 

Here, we propose to invoke the existence of a spatial 
discontinuity of the fine structure constant, and perhaps 
also of other constants, due to the existence of a domain 
wall crossing our Hubble volume. As depicted on Fig. Q] 
(left), if such a domain wall exists the vacuum expec- 
tation value of the scalar <f> that supports the domain 
wall changes sign across this hypersurface. It is then 
clear that if the fine structure constant a is a (non-even) 
function of <fi, then it shall take two values, one a+ (the 
larger value of a) in our neighborhood and the second a_ 
at high redshifts in the direction of the wall; see Fig. [T] 
(right). Indeed, on scales larger than our Hubble volume, 
there exists a stochastic distribution of a taking arbitrar- 
ily one of these two values. Since in such a scenario, <fi is 
expected to have a mass much larger than Hq , it is sta- 



bilized in one of its two vacua so that a has to be strictly 
constant in each patch (this is similar to the landscape 
approach to the cosmological constant problem but on 
scales of order the size of our observable universe) . This 
implies that local constraints [l| on the variation of a 
such as atomic clocks, Oklo and meteoritic dating will be 
trivially satisfied. 
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FIG. 1: (a) A domain wall is assumed to cross our Hubble 
volume. It intersects our past light-cone on a 2-dimensional 
spatial hypersurface characterized by the redshift of the wall 
in a direction n, zw(n)- z, is the lowest redshift at which 
the wall can be observed and corresponds to a direction n*. 
According to Ref. n, should point towards right ascen- 
sion 17.3 ± 0.6 hours and declination —61 ± 9 deg. (b) On a 
constant time hypersurface, the wall cross our Hubble horizon 
so that the fine structure constant takes 2 values: a+ in our 
neighborhood and a_ on the other side of the wall. 

The simplest way to implement this idea is to consider 
the following theory 



S = 



\ M 2 p R-\{d^ 



v&) + \b f {4>)fI u 



^iipjpipj - Bj^mji/jjipj 



-gd 4 x, (3) 



where M~* 
scalar field 
tial 



= 8ttG is the reduced Planck mass. The 
is assumed to have a simple quartic poten- 



1 



2\2 



V{cj>) = -A(^ - rf) 



(4) 



and a coupling to the Faraday tensor of electromagnetism 
as well as to the fermions ipj. One could generalize the 
theory so that (f> couples to other gauge fields as well and 
even to dark matter as considered in Refs. 0, l23l l24j- 
but this would not change our argument. The coupling 
functions B; are assumed to be of the form 



Bi(4>) = exp £ 



1 + 6 



(5) 



1 This could lead to the mechanism now known as chameleon ll 2| . 



where the coefficients £j are constant and M» is a 
mass scale. This model depends on the parameters 
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(A, M*, 77, £f, £i) and we shall assume here that only £p is 
non-vanishing at tree level. Nevertheless, the scalar field 
inevitably couples to nucleons radiatively through £jv = 
rn^(N\^ F /i)F^ v \N)M. This yields £ p - -0.0007£ F 
and £ n = 0.00015£f [25| respectively for the proton and 
neutron. Since most baryons in the universe are protons, 
we shall take £jy = Cp f° r simplicity in our estimates. 

This model is a generalization of that introduced by 
Bekenstein [26] and is useful for the investigation of the 
connection between the cosmological variation of the fun- 
damental constants and the amplitude of the violation 
of free fall in the Solar system 0, [3 H3 and similar 
couplings have been argued to be a generic prediction 
of string theory at low energy ||. The main difference 
between the model studied here and previous models is 
that the scalar field is assumed to be heavy so that it is 
stabilized, hence we do not expect any local violation of 
the equivalence principle. Indeed, the current model does 
not exhibit any temporal variation of constants once the 
phase transition has occurred. 

The evolution of the field is dictated by the Klein- 
Gordon equation (TT]). The effective potential gets three 
main contributions in addition to the potential ((U); 
namely from (i) the coupling of <f> to the electromagnetic 
binding energy of the matter, that is, to pbaryon, from 
{it} loop corrections that will scale as £ F 2 T 4 /M 2 and 
from (m) finite temperature corrections which scale as 
(d 2 V/d<j) 2 )(t) 2 /24: if the field is in equilibrium. Note that 
there is no coupling to the radiation energy density since 
(F 2 ) = 0. Thus, the effective potential has the form 

V eS = V(cf>) + 6v^Pbaryo„ + + (6) 

To determine the typical order of magnitude of the 
model parameters let us first ignore the thermal correc- 
tions and the coupling to matter. To reproduce a change 
in a through the domain wall matching the claimed spa- 
tial variation 11611 . one needs 



Aa n K 
~ 2f F — ~ few x 1(T 6 . 

a M* 



(7) 



For simplicity, we shall assume in our numerical estimates 
that r\ = M„, so that £p ~ (note that £p can be 

chosen positive without loss of generality) . According to 
the claim in Ref. (l6| . we would need to be living in the 
vacuum with greater a, which we denoted a+. A greater 
value of a means a lower value of B Fl which (for positive 
£f) implies that <p = —77 at our location. Since £ p < 0, 
the a+ vacuum has a slightly greater energy density than 
the a,- vacuum. 

Once formed, a static domain wall has a field configura- 
tion <p(z) — 7/tanh(z/2: c ) with z c = (2//j) = (y / \/2rj)~ 1 
being the typical thickness of the wall (this is sim- 
ply the solution of the equation of motion for <fr in 
Minkowski spacetime and with the potential ((H); see 
e.g. Refs. [HI, [29[. It follows that its energy density is 
Pwaii = A77 4 / [2 cosl^lz/zc)] with a surface energy den- 
sity t/waii = 2\/2 Ary 3 /3 obtained by integrating /O wa ii over 



the transverse dimension. For a domain wall spreading 
on a scale H^ 1 , the total energy density is of order UH . 
It follows that the contribution of the wall to the energy 
budget of the universe is of order [30| 



wall 



Uyra.ll Hp 
PO 



V 



100 MeV 



(8) 



where po is the current total energy density of the uni- 
verse and where we have fixed A = 1. In the following 
we shall thus assume 77 = O (MeV), so that the energy 
density in the wall is sufficiently small. It follows that 
the typical values of the parameters of our model are 



A ~ 1, r) ~ 1 MeV, 77 = M,*?, 77 ~ 1, 



and 



£ F f} ~ 10- 



-7 x 10- 4 £ F . 



(9) 



(10) 



Similarly a network of domain walls, which were however 
assumed to be frustrated, with rj ~ 100 keV was con- 
sidered in Ref. [3l| as a possible explanation for the late 
time acceleration of the cosmic expansion. 

Let us now discuss the cosmological evolution associ- 
ated with the potential ([6]) . In writing this potential, we 
have implicitly assumed that the quanta of <fi have the 
same temperature as Standard Model fields. The interac- 
tion rate for the scattering between two photons and two 
quanta of <f> is parametrically given by T ~ £ F T 5 /Mf. 
For our choice ([9lfT0|) of parameters, this process is effec- 
tive (r > H) as long as T > 10 MeV. At lower temper- 
atures the rate of <fi self-interactions remains high, due 
to the much stronger \<j) A vertex. We therefore conclude 
that the potential ((6]) is perfectly justified. 

Ignoring the small corrections proportional to £f and 
£iv, the %2 symmetry is restored at T > Tq = 2r\. The 
third term in (|6|) is then much smaller than the fourth 
at any temperature of interest, and can thus be disre- 
garded. The linear term in <p instead shifts the mini- 
mum to a slightly positive value of <f> during the unbroken 
phase (since £ p < 0), and introduces a tiny discrepancy 
between the potential of the two minima in the broken 
phase. At the phase transition, these effects are negligible 
with respect to the relevant scale, T, of the 8(f) fluctua- 
tions. 2 It is simple to check that the correlation length 
in the standard Kibble-Zurek mechanism [29], £ ~ T^ 1 
is much smaller than the horizon size at the time of the 
transition. The wall thus formed at a typical redshift of 
l + Zi = T c /T where T - 2.348 x 10~ 4 eV is the CMB 



2 For T > Tq , the potential has a minimum at small but positive tf>. 
At T < Tq the potential has two minima, and a maximum. The 
local minimum and the maximum appear at T = Tq, and they 
coincide at that time. Denoting by <f>- their common value at 
Tq , and by rf>+ the value of the true minimum at Tq , we find that 
0+ - 0_ ~ 10- 5 r?(5 F MeV/10- 6 AA/») 1/3 , and that V («£_)- 
V ((p+) ~ 10" 18 V (£ F MeV/10- 6 AM*) 4/3 , where V = A»7 4 /4 . 
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temperature today, so that Zf ~ 8.5 x 10 9 . In particular, 
our model differs from the late time phase transition in 
a proposed in Ref . [HJ , where the characteristic scale in 
the potential is meV. 

The spatial distribution of domain walls that form at 
the phase transition can be studied from percolation the- 
ory (33| , which concludes that, shortly after the transi- 
tion, the system must be dominated by a large wall with 
an extremely complicated structure spreading along the 
entire Universe [34]]. Smaller closed walls are also present, 
but they quickly contract and decay. Typically, the evo- 
lution of the system eventually leads to one large wall 
per Hubble radius. In the case considered here, the two 
minima on different sides of the wall have different en- 
ergy due to the linear term in tf> present in ([5]). As a 
consequence, the wall is subject to a force towards the 
region of higher potential. In our case, this corresponds 
to the wall moving towards our location. 

Due to this, a number of consequences may be expected 
if the wall moves at a relativistic speed today. Firstly, the 
absorption regions, from which the variation of a is de- 
duced, need to be in the a_ vacuum, while we need to be 
in the a+ vacuum. Light passing through these regions 
and arriving to us needs to cross the wall. This is highly 
constrained if the wall itself is moving towards us at rel- 
ativistic speed. Secondly, photons crossing the wall can 
be reflected by it with an O (8a/ a) 2 probability (this can 
be easily obtained by computing the reflection and trans- 
mission coefficients of a flux of photons incident on the 
wall); a further (£tv/£f) 2 suppression is present for mat- 
ter. Even if it is a small probability, the reflected photon 
(or nucleon) will have a large energy if the 7 factor of the 
wall is high. This could lead to phenomenological signa- 
tures when, for example, the wall crosses a star. Thirdly, 
the angle of the cone in which the cone is visible, Eq. (TIT)) 
below, is affected if the wall moves relativistically. Note, 
however, that interestingly the term in Bp(<fr)F 2 does 
not affect the equation of propagation of photons in the 
eikonal approximation at first order [351 ] . 

The motion of the wall can be derived from the dynam- 
ics of extended objects [36| according to T^^K^ = f 
where an( i are the stress-energy tensor and ex- 
trinsic curvature of the wall and / the force acting on it. 
This reduces to 



d 2 x 3 
~dr r 



dx^ dx v 7 AV 



dr dr R U wa n 



(11) 



for the motion of a flat wall in the transverse direction 
x 3 = z, where r is the proper time measured by an ob- 
server on the wall, R the scale factor of the universe, 
= 2?7£jvpbaryon/Af* is the potential difference be- 
tween the two sides (see Ref. [37| when the effects of 
wall curvature are relevant). As x 3 is a comoving co- 
ordinate, the physical velocity of the wall is given by 
v = Rdx 3 /dt = (R/j) dx 3 /dr where t is the physical 
time, and 7 = dt/dr = l/yl — 



quantities, Eq. (|11[) can be rewritten as 



(12) 



From the expressions for AV and t/ wa ii, and from the 
current baryon density Pbaryon.o — 1.8 x 10~ 48 GeV , we 
deduce that F — 2.7 x 1CT 48 (1 + z f f GeV with 



10" 



(_5_y 

VI MeV J 



/A- 1 / 2 . 



(13) 



This equation of motion can be solved numerically as- 
suming that the wall starts at rest at T = Tq = 2 MeV, 3 
and for / = 1; the result is depicted in Fig. [2] We see 
that, right after its formation, the wall accelerates to a 
large boost factor, 7^0 (10 6 ). However, it gradually 
slows down due to Hubble friction. The current peculiar 
velocity of the wall is vq ~ 0.004, so that none of the 
effects mentioned above is an issue. We also see that, in 
the non-relativistic regime, the velocity of the wall scales 
linearly with the parameter /. 




10000 



1e+06 
z+1 



1e+08 



1e+10 



FIG. 2: Evolution of the physical velocity of the wall v (solid 
line) and of the associated boost factor 7 (dotted line), plotted 
as a function of cosmological redshift z. We recall that time 
evolves from right to left in the Figure. 

To be viable, our model must satisfy an additional set 
of constraints that we now summarize. 

1. CMB constraints. While the effects of cosmic 
strings on the CMB have been extensively stud- 
ied [381 ] , domain walls have not been widely consid- 
ered since they were thought to be formed at much 
higher energy and thus have a dramatic effect on 
the CMB. 



In terms of physical 



3 We can show analytically that the velocity of the wall today 
depends only logarithmically on the initial temperature. 
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A single wall would contribute to the tempera- 
ture anisotropy via the integrated Sachs Wolfe ef- 
fect [1^]. For a static universe, and in the non- 
rclativistic regime, the gravitational potential gen- 
erated by a very large (planar) wall grows linearly 
with the distance L from the wall (notice that this 
has a zero net effect if the distance between the 
source and the wall is equal to that between the 
wall and the observer). For a typical cosmological 
distance ~ Hq , we estimate ST/T ~ G U wa u H^ 1 . 
This results in 



J / CMB 



10" 



VI MeV/ 



(14) 



which constrains r\ to be smaller than a few MeV. 
Another source of temperature anisotropy is related 
to the fact that a is not constant across the visi- 
ble universe. This is not an issue for the model 
we are discussing, since the cmb data are only sen- 
sitive to a variation Aa/a ~ O (10~ 2 ) or greater 
40]. Finally, we can also neglect the temperature 
anisotropy related to the probability that a CMB 
photon is not transmitted through the wall, which 
is of O (Aa/a) 2 , as we have already mentioned. 

2. Astrophysical constraints. Although our scalar is 
relatively heavy, ~ 1 MeV, it can be produced 
in supernovae. The production rate of scalars 
through inverse decay is roughly, 



- 77->0 1 



Ml ' 



(15) 



In principle for scalars with mass < T, this 
could result in an excessive energy loss rate. How- 
ever, these scalars decay to two photons with a rate 
~ £f.///M* ~ £f*M*. Requiring that their de- 
cay length is smaller than the size of the core leads 
to a lower bound on their mass 



AL 



lMeV 



> 0(10" 



10- 



(16) 



(for typical energies of order T ~ 30 MeV). Thus, 
for our choice ([9l fTU|) of parameters, these scalars 
decay within the core and there is no energy loss. 

3. Tunnelling to the true vacuum. Contrary to a stan- 
dard domain wall, the non-minimal coupling in- 
duces a shift between the two minima. The lifetime 
of the false vacuum is of the order 141] 



Aexp 



27tt 2 gg 
8 AV 3 



with S = fl ^fV{(j))<l(j) ~ 2t7v / AV. Now, with 
AV ~ ^N^moPcntiv /M*)(l+z) 3 and assuming A — 
M" 1 , we conclude that r = Hq 1 f T with 



f H ° 

U = MZ^ 



54tt 2 M 4 

Qmo£,N Pcrit 



r) 3 (l + z)- 



The argument in the exponential is always very 
large and scales as 4.5 x 10 54 ?) 3 (1 + z)~ 3 . Even at 
the time the wall is formed, Zf ~ 10 10 , the factor is 
larger than 10 24 . This means that the wall forms at 
a time where the false vacuum has a lifetime larger 
than our Hubble time. Since the lifetime of the wall 
increases with time, we are guaranteed that today, 
the wall is effectively stable. 

Our model also has some specific observational predic- 
tions that arise from the fact that a takes two discrete 
values. Let us denote n* the direction of the wall and 
z* the redshift of its closest position and x* = x( z *) the 
comoving radial distance to which it corresponds. Since 
the equation of our past-light cone is x — Vo — V an d that 
of the wall x — X*l cos# if it is assumed to be at rest in 
the cosmological rest frame. It follows that the redshift of 
the wall in a direction n is given by x( z ) = X*l cos# with 
cos = n* n. Solving this equation gives the dependence 
of a with z and 6 and is plotted in Fig. [3l Moreover since 
our observable universe as a finite radius, the discontinu- 
ity can be observed only in a cone of angle 81 around the 
direction n* with 



cos 



Qi = x*/xn- 



(17) 



61 depends on z* 
Typically, cos 6*/ 



and on the cosmological parameters. 
0.345 if z* - 1.8. 




FIG. 3: Summary of the observable prediction of our model. 
Assuming the wall is at a redshift z* ~ 1.8 at its closest posi- 
tion to us in a direction n„, then the fine structure constant 
is a function of cos 6 — n» • n and z. The blue region corre- 
sponds to a — o_ and the upper region to a = a+. The fact 
that our observable universe has a finite radius implies that 
we shall detect no variation for angle larger than 81. 

In summary, we have proposed a two vacuum solution 
to produce a large scale spatial variation of the fine struc- 
ture constant, which is not associated to a local time 
variation. In the model presented, the scalar field was 
coupled to electromagnetism and to nucleons through 
the photon contribution to the nucleon mass. It is of 
course quite plausible that this coupling extends to other 
gauge fields (and perhaps Yukawa couplings) as well. In 
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that case, we should expect not only spatial variations in 
a, but also coupled spatial variations in other quantities 
such as particle masses and Aqcd [HI 53 ■ 

For example, one might expect a variation in the 
proton-to-electron mass ratio, /i. If due to coupled vari- 
ations of fundamental constants, one typically expects 
that it is correlated to the variation of a; A/i//i ~ 
— 50Aa/a [431 ] - Molecular hydrogen transitions in quasar 
absor ptio n systems yield a limit A/i//i = (—2.6 ± 3.0) x 
10 -6 [44j. Searches for spatial variations in the proton 
to electron mass ratio in the Milky Way produced a 
non-zero result of A/i/n = (2.2 ± 0.4 ± 0.3) x 10" 8 and 
conservatively implies an upper limit of 3 x 10~ 8 . Analo- 
gous searches for a spatial variation of a 2 // yield an upper 
limit of 3.7 x 45]. Taken together, these would im- 

ply an upper limit of Aa/a < 2 x 10~ 7 . In the context 
of the model presented here, this limit would present no 
difficulty with the Keck/Hires observations indicating a 
variation in a. The Milky Way being entirely in the a+ 
vacuum would show no variations in any of these quanti- 
ties and these observations therefore could not constrain 
variations in the a_ vacuum. 

We also note that we might expect a small effect on 
the light element abundances produced in big bang nucle- 
osynthesis [40]. Whether or not the transition happened 
before or after BBN, the average light element abundance 
should correspond to abundances using constants deter- 
mined at 4> = (or at a slightly smaller value of a). How- 



ever variations at the level of 10 -6 would hardly be per- 
ceptible in element abundances. Similarly a small change 
in a would very slightly affect the CMB through recom- 
bination but this too would be imperceptible. The 
same may be true for future constraints that can be set 
from signals originating from the absorption of the CMB 
at 21 cm hyperfine transition of the neutral atomic hy- 
drogen [47j or from stellar evolution [48| with sensitivities 
which are typically of the order of 10 -2 for z > 30 and 
10 -5 for z ~ 15 — 20 respectively. 

Other standard powerful systems [l[ for testing the 
variation in a are the Oklo phenomenon and meteoritic 
studies. But as these are both on our side of the wall, we 
would expect that they are too local to probe any vari- 
ation in a. Similarly, we would not expect any positive 
result from atomic clock measurements of changes in a 
unless of course the wall passes through the solar system 
in the near future. Given the sensitivity of these exper- 
iments (see e.g., Ref. [Hj]), there would be no mistaking 
the event. 
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